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ABSTRACT 


The state variable formulation of a Foster Wheeler 
ESD-[boiler is developed from fundamental principles. 
The response of the model for various input signals is 
determined using CSMP -],the IBM simulation language. 
The sensitivity of the model to various coefficient 
values is noted as are the characteristics of various 


system states for small perturbation values. 
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NOMENCLATURE 


As, Ay,An = Superheater, riser and downcomer cross-sectional 
2 
areas respectively (ft ) 
C = evaporation ievel proportionality constant 
(rt-sec/Lb) 
C.,Cy,Cp = heat capacitance for superheater, riser and downcomer 


tubes respectively (BTU/1b-8) 
Cg = average heat capacitance of combustion gases tines 
air-fuel ratio (BTU/1b-R) 
Ci = heat capacitance for feedwater (BTU/1b-2) 
Cys= average heat capacitance of superneated steam 
(BTU /1b-R) 
§.) fy, fp = friction coefficients for Superheater, riser and 


2 5 
downcomer tubes, respectively (sec /ft ) 


2 
= acceleration due to gravity (ft/sec ) 


a= entnalpy of saturated vapour corresponding to P 


(BIU/1b) 
hue= enthalpy of saturated liguid corresponding to - 


(BTU/1b) 
hye= enthalpy of liguii in mud-drunm (BTU/1b) 
hy = enthalpy of drum and downcomer liguid (BTU/1b) 
hy = enthalpy of evaporation corresponding to Ea (BTU/1b) 


K,,= air cooling coefficient at superheater bank (BTU/sec) 
Kay= air cooling coefficient at riser bank (BTU/sec) 
Ky. Kg = heat~-transfer coefficients from combustion gas to 
Superheater tubes, and from superheater tubes 
to steam, respectively (BTU/1b-@) 


Kga,Ky) Ky = heat-transter coefficients from combustion gas to 








Ni 


cviser tubes and from riser tubes to boiling 


® ° 6 
liguid respectively (BTU/1b-R), (BTU/1b-R ), (BIU/1b-2) 
Evaporation rate constant of drum liguid (Lb/sec-R) 


constants ror state equations of saturated steam 


2 =1 
(keto 7lb), (ft) 


Superheater, riser and downcomer tube lengths 
respectively (ft) 

mass of drum liguid (1b) 

mass of Superheater tubes (1b) 


mass of riser tubes (15) 


2 
drum pressure (lb/ft ) 


2 
Superheater outlet pressure (l1byft ) 


2 
mud-drum pressure (lb/ft ) 


heat-input rate from tube walls into the superheated 
steam (BIU/sec) 


= heat-input rate from hot gasses into superheater tube 


Walls (BTU/sec) 

heat-input rate from riser tube walls into boiling 
liguid (BTU/sec) 

heat-input rate from hot gasses into riser tube walls 
(BTU/sec) 

reedwater temperature (R) 


Saturation temperature corresponding to P (R) 
B 


Superheater outlet temperature (R) 
drum and downcomer liquid temperature (R) 
Superheater tube-wall temperature (R) 
riser tube-wall temperature (RB) 
average gaS temperature at Superheater banks (R) 
average gas temperature at riser banks (R) 
compustion gasS temperature entering sSuperheater banks 
(R) 
3 
volume of vapor phase in drum (ft ) 


10 








Vy 


= saturated vapor density corresponding to ae 


velocity of riser mixture (ft/sec) 


= velocity of downcomer water (ft/sec) 


3 
total drum volume (ft ) 


steam mass-flow rate at the superheater outlet 
(lbE/sec) 
fuel mass-flow rate (1b/sec) 


= feedwater mass-flow rate (lb/sec) 


downcomer mass-flow rate (1bd/sec) 


Ciser mass-flow rate (l1b/sec) 


= air mass-flow rate from blower (1D/sec) 


= chemically correct +50% excess air rate (l1b/sec) 


mass-evaporation rate from drum liguid surface 
(lb/sec) 

steam mass-flow rate from drum into superheater 
(Lb/seéec) 

guality of mixture leaving riser 

throttle opening (4%) 


drum liguid level (ft) 


3 
Sup2rheater outlet density (lb/ft ) 
saturated liguid density corresponding to P 
density of liquid vapor mixture leaving riser 


3 
(lb/ft ) 


evaporation level (ft) 
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3 
(lb/ft ) 
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I. INTRODUCTION 


Boilers as understood by marine engineers are closed 
vessels containing water which by the application of neat is 
converted into steam at any designed pressure. This steam 
is then used for the production, through machinery, of 
useful work. A dynamic model of steam turbine machinery 
consists of a boiler model and a turbine model. Tne 
difficult part is the boiler, as a load change causes 
Variations in somé important properties such as koiler 
pressure, temperature and drum water level.The analysis in 
this paper will use Chiens dynamic analysis (1) of a boiler 
as a reference.Chien considered a naval boiler which for 
purposes of analysis was divided into four sections namely a 


Superheater, a downcomer-riser loop, a drum anda gas pata. 


The principles of thermodynamics,heat-transfer and fluid 
mechanics were used to describe the dynamic behaviour 


corresponding to each section of the boiler and these were 


derived from equations of continuity, energy, 
heat-transfer,and momentun. The eguations involve partial 
differentials as well as nonlinear terms. These equations 


were reduced to the ordinary linear equation form by 
applying small perturbation and difference equation 
technigues.Linear equations thus obtained were reduced to 
ten state variable equations and solved by digital computer 


techniques. 


Since there iS an increaSing interest in toiler 
modelling, the objective of this thesis was to develop a 
comprehensive boiler Simulation modeit in a form useful for 


modern control (i.¢e., multivariable control) analysis. 
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II. BOILER MODsL 


A. BOILER CONSIDERATIONS 


The control problems of high pressure boilers have 
become more and more critical, both from the operational and 
the economical points of view. A dynamic analysis is the 
method to be used for a control-system analysis. ie 
consists of a complete understanding of the process to be 
controlled and the effects of physical and chemical changes. 
The analysis is not exact by any means but the results 
obtained should be in good qualitative agreement with actual 
tests. The major difficulty in boiler analysis is the fact 
that the whole system is very complex and contains numerous 


Variables. 


Chien(1) considered a naval boiler which for purposes of 
analysis was divided into four distinct sections namely a 
superheater, a downcomer-riser loop, a drum anda gas path. 
The detail of analysis is described in Part B of this 
section. The basic equations used jn the analysis are those 
of continuity,energy (heat-transfer),momentum and the state 
equations. These eguations involve partial differentiation 
as well as nonlinearities. Generally,the equations have tae 


form 


4 (x i ae) an) 200 


To eliminate the nonlinearities one uses perturbation 
theory, which effectively approximates the response of the 


4 
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system to small signal changes about a chosen operating 


Sonditi1on. Thus, the equation is perturbed about its steady 


state operating condition to give the linearized forn. 


Hence it can be written as 


Of ax, May, : 
+ AZ+---- = 0 
ox vy Y Oz 
The perturbed variables are AX, AY, AZ, etc. and the 
partial differentials ( Sf, Of, Af, etc.) that form the 
aH x Oy Oz 
coefficients of the perturbed variables are evaluated at 


steady state operating conditions. This technique was also 
followed by Whalley (2) in modeling the same boiler plant. 
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B. BOILER ANALYSIS 


The boiler studied in this analysis is a Foster Wheeler 
D-type marine boiler; which is an oil fired,two-drum,natural 
circulation unit having a rated output of 28,800 lbs/hour at 
350 iss 78 Cameo. 1200 F. A cross sectional arrangement of 
the unit is shown in Figure 1.The following assumptions 
apply to the physical simplifications in each of four 


sections. 


= =a = ae Se SS ees SS 


(a) The inertial effects of superheated steam are 
neglected. 

(b) The superheater tubes are assumed to be a single 
capacitance with restriction on the drum side and another 
restriction on the load side. 


(c) DesSuperheaters are not considered. 


(a) Only natural circulation exists. 

(b) No boiling takes place in the downcomers. 

(c) Vapor and tLliguid velocities in the riser are 
identical. 

(dj Heat=transter tates to the boiling liquid from the 
tube walls are proportional to the cubve of temperature 
difference between the wall and the liguid. 

(e) Steam quality is uniform in the riser. 

(tf) Liguid temperature is always the same as the 
Saturation temperature corresponding to drum pressure. 

(g) Downcomer liguid temperature is the same as tne 


drum liyguid temperature. 
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3. Drum 


(a) There iS notemperature gradient across the drun 
vapor phase,and the temperature is always the saturation 
temperature corresponding to drum pressure. 

(b) The liguid phase has no temperature gradient other 
than acrosS a very thin boundary layer at the drum surface. 

. (c) Evaporation or condensation rate in the drum is 
proportional to the difference or liguid and saturation 
temperatures. 

(d) Feedwater temperature is assumed to be constant. 

(e) Liquid-level changes due to bubble formation in the 


drum are neglected. 


4. Gas Path 





(a) The air~fuel ratio is assumed to be constant. 

(b) Temperature of combustion gas entering superheater 
is proportional to the firing rate. 

(c) Waterwalls are lumped with the riser-banks. 

(d) The heat-transfer rate at e?ach tube bank is 
determined by the tube wall temperature and the average gas 
temperature. 

(e) Inertia of the hot gases is neglected, that is, 
velocity changes take place instantaneously. 

(£) Delays due to the heat capacitance of the hot gases 
are neglected, that 1s, temperature changes take place 
instantaneously in combustion gases. 

(g) All heat transfer is due to turbulent convection 


and radiation. 


Th2 following steps are taken in developing the 
equations for each of the four sections. A simple Schematic 
diagram of each sub-section is included at the beginning of 
each part of the analysis. A brief statement of the physical 
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Situation is included under the headings. The dummy 


coefficients such as aus Die Ci Dat were used for 


convenience),and their values are declared in the Appendices. 


The quantities such as wae Bae nee etc. are the steady state 


values of W, hh, Ba respectively. These values can be 
B Ss 


evaluated trom the original unperturbed nonlinear equations 


by setting all derivative terms to zero and solving for the 
unknown values in term of known values. 
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Figure 2 - SCHEMATIC DIAGRAM OF SUPERHEATER 


ify, 





Mass flow in - Mass flow out = Rate of change of enclosed 


nass 


Ne Ws 


G (bss Pa) 


The perturbed eguation is: 


OL; 


O, AW, +4 AW, = OSAP, (1) 


where S§ is the' Laplace variable denoting differentiation 


With respect to time 


Foal 


Oo. Superheater energy balance 


Heat input from tube wall + Heat in incoming steam - Heat 


With outgoing steam = Rate of change of internal energy 
d 
Q.+W,h, WAL = (shes) 
The perturbed eguation is ; 
Aa +Wah, +h Ws Oh, hw, = LA. (P.SAh+ h SAP, | 


Since, Ah,=0, Ah.=C, ATs and LAsS AP, = a AW 


8 S 
After substitution and rearrangement: 
(LAs PeCns} SAT, us (h,-h, )AW, ~ AQ, = “WoC yg AT, (2) 


OL, 


0,5 AT, + @,AW,+0,AQ, = GAT, 
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Heat input to wall - Heat output from wall = Rate of change of 
energy stored in superheater wall 


GS S SW 


d 
Q = = at Nest ) 


The perturbed equation is: 


AQgs — AG, MoCoS ATS 


OL, 


e. 
> 
a) 
+ 
2 
ae 
2 
il 
2 
W 
> 
—, 


3 4s g S 10 SW (3 ) 


The following empirical formula includes the cooling effect 


if air supply W,differs from ideal air supply We (Whalley (2)): 


0.6 2 
Qos = Kgs Wy (Tys - Tow) - Kas ( _ Wa ) 


The perturbed equation is: 


0.6 Ky, oF (Wac7 Wa 
Aa, = Sah (Tae “7 JAW, # Ky Tig (ATye= A Tgy) 2K, eH) a 


a 
Wace 


Or; 


sae = Tie AWE + O,,CATI.- AT. ) +0, AW, (4) 
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a ee a ae ee ee es ee ee ee =e So ES ee ee a > oe 


Tae perturbed equation is: 
0.3 = = — 0.8 
§ ee 3 \ ‘sw T,) AW, + ac (AT. - AT, ) 


OL, 


L5 Se Og AM, us ig (OTe = Ate) (5) 


10. Heat transfer from combustion gas 


= 2 SSS SS ee eS ee a eS ey See ae ES Se ES 


Average gas temperature at sSuperheater wall is T where 
gs 
0.5 Qs, 
ea" 
and , is related to the heat capacitance of the combustion 
gas by : 


Tos = Te - 


where cis ccnstant 
Thus: 


Qo. = 2C.(Wy + Bion) @teaa) 


The perturbed equation is; 


Aa. - 2c Wi + Wi) AT. + (T.-T ye) OW, + (Te ~Tys AW, (Wy + Wa) A734 


22 





Rearranging: 


+ AG 


2c. a6 (We +W,) AT. +(T. ~Tyg) AWg + (Te -Ty.) AW = (Wp +) Tr, 


OL, 





Since : 





We = fa As Va 
Z 
L 
B 
$7 905 | Paagl /® 
Substitute: f{ _ bls 
a 
The perturbed equation is 
- ee 
Ap -~ AP, — 7 Wa AW We A 
B 35 $= a -ft. = 
Pe Nas Ie 


Rearranging : 


. 


= 2 
Wa A: a 
is oe AMe + AP “he AP, = APs 


8 
OL, 


a, AY, + 0,, AP, + a = @, AR, 
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Q = 
8 4s Big AT; + on AWS - a, 40, + 1241, 


(6) 


7) 
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Figure 3 - SCHEMATIC DIAGRAM OF RISER-DOWNCOMG&R LOOP 





Mass flow into riser - Mass flow out = Rate of change of 


enclosed mass 
WW _ 34 ols coe 1) (+P, )/2| 


The perturbed equation is: 


AW, ~ AW = L] Ay Ly SAP ~Ay(P+Py) SOY + Ar SAP| 


since 
de fs Pw 

then, 2 
A oe ax+L%XA 
Pf PUR, zs) 5 = 


If the starting point of evaporation in the riser tube is y 


above the water drum exit where ;: 


Wo = {-CWy 


then, 
Ay = Se Sale 


wer y S Ie ‘can be neglected since » 1s small 
Substituting (a and - and rearranging : 


o.5A hep (4, J, ) SAX 40-5 Ay LyX K. (4. sar, 
FOS (P+ Py)CSAW, = Aw,~ AW 


or, 


— 


GQ, SAX+G, SAR +A, SAW, = A,,AWy + O35, AW (8) 
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13. Riser momentum equation 


Pressure drop - frictional losses - gravity head - 
exit loss - monentum flux change = Inertia force 


4 
fe 


fp by W. W 
Par Seen = ALP = 2... DA 
Ue P, ) Ay. Ar gdp yr ¥ ¥ vp 24 (pray? / 


9 
on 


-| stems ae ~ te du 
J (P*y) 9 (Py Av)? 4 dt 
Thereiore: 


a) LCS 5 + Ba - gee 


After perturkation and substitution for Ap as done before 





in the Riser Continuity Equation, the perturbed equation is: 








Vea Pw Pw B W 
FX 7 ; 
t 23 ae 4 ie ($-¢)fay = =" S AW 
Ps SPs Ss fa Pu Pe qv 
where (Set . : 
Dy 2 
—" , 
Ay 9 
Rearranging - 





2Ww 
AR, OP, - (b,+0,8 AW+ by AX- =F AW, #o< OA, 
Vi] OW 
where 








See _X Hi (in, 8 yal and is small enough 
Ps Aw\ Dye 2 fa u to be neglected 


Hire LeLOne ¢ 


AP,- AP, = AW+ a SAW+b, Ax- 2 AW, 
TP 


OL, 


QAP, +O,, AP, = 0,,AW4 0, ,SAW+0, AX +0, AW (3) 


Heat From hot gas to riser tubes-Heat from tubes to mixture 


=Rate of change of internal energy stored in tubes 


d , 
Sas = Q = dt CRE) 


The perturbed equation is ; 


AQ,. SG swe So tla a 


OL, 


0 ee 
yee By AT ay (10) 
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15. Riser heat balance equation 


Heat input to risers-Heat output from risers 


=Rate of change of internal energy 


Q.- Wyhyy- WA = ¢. | 0.5 Ar(te-y) (P49, Vh | 


From the continuity equation ; 
d 5 is Pa — —~ 
go Anlte CP + Pa) = y= 
After substitution and perturbation : 


AQ + WAR + Ch, -n) Aw, = Ah =- 0.5 Ay ew (P+ P,) SAh 


Since: 


bh = Me, + Ng 
If changes in latent heat Ahi are neglected,then 
Ah = Ahn + hea AX 


From steam tables : 


Ah, = KAP 


A hens = At, 
Therefore; 


0-5 Arlo (P+ py) KeS AP, 40-5 Anly (P+ Py) hg SAX ~ Aa 
= (hus ~ h)awy - WK ~ Whe Ax +h, AT, 


Or, 


q P +a : 
SAP, +O,,SAX4+0, Aa. =a, Aw + ART AX+0 At, i) 
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16. Riser heat transfer equation 


An empirical formula that accounts for turbulent gas flow 


is adopted and gives: 


The perturbed equation is: 


0.6 Kan (Ton - Tr = 
A Las = ae — AWS 4 Kap Wt CATs3—AT,.,) 
NE 


» 2Kaw(Wac = Wa) AW, + 4ky (T. 


ie Wat 


ga AT a Ey, ») 
Rearranging . 
8 6.6 a 
= a 0, : = ) 
AG... Ngo Nf AT gs = — Kgs Wt ATay + ei Tas Taw AWs 
We 
t 


2 Ko (Wac-W ) _ 3 
wo ae AW, + 4K, 79, ATos 
AC 


~ 


hence ° 


aK, sii “+ AKT Taq) AT, qn =~ (Kost + 4k) F, has 


ee W%6Kg, (Tga -Taw) A We 


Tito 


+ 2Kaw(Woe - Wa) AW, 
Wor 


Or, 


A yg AQga + Teg ATa, = AAT ey +O, AM +O. AW, Cle) 
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Lue iser heat transfer to boliing fluid 


The heat transfer rate is assumed to be proportional to the 
cube of temperature difference between wall and mixture. 


Continuing to follow Whalley (2): 


OQ, > - (oa 
The perturbed equation is 
e 2 
AQ = 3 Ka (Tay- Ta) (AT. ~ ATS) 
CL, 
a5, 48, = Os ATaw + O55 ATy Cle 


18. Downcomer momentum equation 


>» Forces due to pressure difference~frictional loss 
+gravitational head-entrance loss = Inertia force 


Hence ¢ 





_ apy on ne Ma) 
an iy = a, (Va) Ay Ao + Pw Ap (bp+2) — 7g Puts 


( LyAy f, Vw 





Since; Ww = Ay Ap Vn 
(P,- Pa) pty Wu Wi a boy 
7 9 fis Aa Pave 29 AR dt QAn * 
2 
( P, ~ Py) (e+ os) ! | qe Oe 4 (Wy) 
D oP 2 ia dt 


30 





The perturbed equation is; 


AP, - AR, = (by+ay,S)AW, —L, AA, 





wheres: . 
byw = zie ( toss Hp 40.5) 
GPudp \ Po 
Ow = mee) 
Q Ab 


Since flow is incompressible Can constant both in 
and space),then AP a0 
Therefore : 


19. Drum steam mass balance equation 


> Vapor rates to/from drum = Rate of change of vapor 
in drun 
Thus , 
Evaporation rate + vapor rate from risers-stean 
d 
rate from drum = “ (Vp 
me) 
Therefore: 
We+XW-W _ 4 ) 
= , = B/a 
ae 


The perturbed equation is; 
We +X WAX — _y 5 
AWse +X AW+WAX AW, = VeS AA +P, SAV, 


Substituting for :! 


A Va = ~AAdy 
A Pa = Ky AP, 
A We = Ke (AT, =i) 
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time 


(14) 


hass 





gives ; 


Yaka SAR - PASAY « Ke CAT, ~AT.) +X Aw+WAX- AW, 


OL, 


Beg S APL + Agi SAY= Aca (ATy-AT,) +O, AW+, AX +A AW, CS) 


%. Heat to and from drum with liquid = Rate of change of the 


2p 
internal energy of the liguid 

Hence: 

Internal energy of riser liguid + heat transfered in 
feedwater - heat taken in by downcomer Liguid - latent heat 


released with evaporating liquid = d (drum liquid internal 
energy) dt 
Therefore: 


Cet eex WT. + Cy WiT; —CyWyTy — Wel, = & (cont) 
But Cy Ci , Cy 2d. 
The perturbed equation is : 
(1-%)Ta AW+ (4~K) WAT, ~ WT, AX 4 WAT, +7; AW, — Wy AT, 
—Ty AWy- We dh, —~ ha A We = MSAT, +T, SAM 


Substitutiag for: 
AWe = Ke (AT, ~ AT, ) 
A lilesee Ko AP, 
Lo we APyAy 


and ie constant,so AT, s 0 


Bi 





Thus; 
(14-X)T, AW4 (1-X) WAT, — WT, AX a 4) WwATy — Ty AWy -Weke APs 


Now Ke (AT,~ AT,) = MSATY HAP, Ty SAY 
Ol, 


Keg OW + Ag7 ATg + eg AX t+ Meg AW; + Ar ATy +07, AW +0, AP, 


+4,,(AW-ATR) = O,SAT, + SAY (416) 


21. Drum liguid mass balance equation 


>, Liguid input/output rates = Rate of change enclosed mass 
of liguid 
Hence; 
feedwater input rate + riser liquid input rate-downcomer 
liguid rate-evaporation rate = am) 

C 


Therefore 


Wit (L-xX) we, - We = Om) 


The perturbed equation is: 


AW; +(1-xX)AwW-WOX- AW. - AW, = SAM 
Substituting for: 
AWe = Ke (ATy- ATa) 


AM = APy Ay 


gives: 


-Ke AT, +A PwSAy 


(1-K) AW-WAX-AWy- Ke AT, + AW; 
OV, 


Ag, AT, + a. SAY 


Qa. AW (7 
83 + QO, 4Xt OA Mt GAT, +o AW, ey) 


33 





22. Equation from throttle valve 


We es Hil Sl )T, ) 
The perturbed equation is : 


AW. = 3 a xy » of AP, + 95; AT, 


Pe OTs 
Ol, 
Og, 45 = AgsAXy + O46 +0, AT, (18) 
23. State equations 
At, = 27 AP K_ AP 
B SB & Bo a B 
B 
OT, 
: 19 
AT. ” a4, AP, 
Ap. = Leap K, AP 
as “ERP 
B 
OL, 
APs = %, AP, os 
? : 
A Tys = AT gs oad 


dhe am Kea 


AN gy = 
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OL, 





A Maw = Qo Nae (22) 
Pe = ih eae We) 
A Pe = cis AT, + <= ofc APs 
dfs 
OL, 
A Pe sy Qi5y ATs +10, APS (23) 
ATco = Kp AW 
Or, 
hale fe ae AW (24) 
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C. STATE VARIABLE EQUATIONS 


The following procedure is an algebraic method to reduce 
the twenty-four governing equations to ten state variable 
equations,that is, in standard matrix form of X (t) = AX(t) + 
BU (t) 

Here X(t) is the matrix of system state variables and 


U(t) is the matrix of system inputs. 


From (1) 
a. Sas = a, AW, a a, AM, 
5 . Saw 4 22 Aw 
APs ms 1 A2 S 
OF» 
AR, > b AW, + b, AW, (1.4) 
From (2) 
a, ATs 27 G,8W,—-a,Aa + a AT. 
“4 Os AW Ce A@ Q7 he 
Ty. =-— — = $+ — Al 1.2 
AT, --@ . ( 
From (5) 
- Se Aw, + Sz AT — 247 AT (1.3) 
AQ. ae ) ae SW we Ss 


Substituting © from (1.3) in (1.2) and rearranging ; 
Ss 


é OQig: 4 a » , Ae. a 
aye = peat aU + pal AT. aes Os oe Ais A16 ae > L7 At, 


U4 Qis Os O4 Aa Ais 41Oy5 SW 
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Ol, 





= 34 
From (3) 
QO te ee = a,44,. + QO. AQ, 
AT, - 28AQ 498 AQ 
=e Gig qs wane 5 
From (4) 
AQ = S28 Aw, + 2 AT, ~ S3AT 4, ay 
q5 Gas sy eed ee” . 
From (6) 
_ O19 AT. 4 220 Aw, + S2t Aw, + 222 AT. 


Setting (1.6) = (1.7) and rearranging: 


ra a2 Ses 
iciewig ss ee aT 


warn Gp Gas 
Qi, G13 Qi; Ais 


low 


Qr9 


(0,,- | AW. + i AT. 


es 
ig _ Azz G13 _ Gaz 
Qi: Oia Ati ayy 


Substituting AT. = Q, Awe fron (24) into (1.8) 
o2 
rearranging: 


ar 





oo Gli > obs Qa 


= re | Qn, 
Q dq, a Q ois teas 
ee ee AN ey aL 
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(4.4) 


iss} 


Ci. 6) 


(1.7) 


(1.8) 


and 





Ol, 


AT; = b AW + bat. —b, AW, : (1.9) 


Substituting ATqs from (1.9) into (1.6) and rearranging : 


: : Q , 
AQ = es + sat) Ang (8s ) AT. + Pa P80 Jay, 


95 \Qyy O14 ry Gis 14 


Or, 


AQ; = b, AW + b AT +b AW, | ( 1. 10) 


Substituting 44,. from (1.10) and AQ, from (1.3) into 


(1.5) and rearranging ; 


AT. seco AWS + eran 2] AT. +(e AW, 


O40 Oto Digi Ags Qs 





AN fhe | AT 
Gis: “1s ° A190 Sys : 
GE, 
AT., = b AN : b. AT Si 2 bg AN, + by. AW, = De AT. (1. it) 
From (17) 


9 a a 
Ay «-—* AT, + 3 Aw + 283 Ax 


O77 77 Az77 
Oo 0, : 

» 2 AW, + —2t At, 4 [32 AW; (4.42) 
177 X77 a77 
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Substituting AT, from (19) into (1.12) 








° a 
Ay = —- 272 Ap, 4 88 Awy S88 AX 
aAr77 : 277 77 
a77 a77 Az; 
Or, 
AY 2-b Ar + by Awe bg AX4 dD, AW, +b, ATy + b,, Aw, fad 
From (16) 
ATy = — 27% AY S68 Aws (677%) at, 4 S68 ax 
Ar, O74 Arq A74 
Or Aw 4 [72 ap, (1.14) 


Q_+4 bes 
at 6s AW: + (a 73) AT, + 


W 
074 O74 O74 


O74 
> 
Substituting Ay from (1.13) and AT, from (19) into (1.14) 


and rearranging : 





; 
Chg as 75, Diy 4 A97-%G7 Igy A735 AR, a SOE Ne O75, Dig AW 
Ar4 Qr4 Az, Ary O74 
a Q : b a Ass, Se 
Peg Bo eg eee eA Cy <2 - 82) Ay 
a. + 0 Q + a Q ES 
4 | 22 eh cl et AT “69 75: 930 AW: 
Ary a7, Ar74 Ar, 
OL, 


AT, = b AP “ b,, AN +b, AX + b, AW, rl b AT, +b, AW; (i,15) 
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From (15) 
AP ee a AY 2 AT, < Ot Ar, 4 S82 Awe Set Axs 6s Aw Clee) 
. O &4 O60 O60 Q 


60 Qéo Ogg 8 
Substituting AY from (1.13) into (1.16) and rearranging - 


Qe. Qe. a a Qeieb a 1a ae 
ap, « {2¢irbsr_ Seater) np 4 (Bea _ Sorbus ay | (M6 _ Ser-bis) 
a 69 A 60 6 460 2 6 Qe 


6a 
_ Ogi» boo AN, fe Qos i Qi. bay Aa, 
D bo Q bo 
We iach QA, ~ 
eA Wieceee ea AW 
Q bo Abo & 


AND SCARE en b, AW + b., AX - b,, dw, 








OT, 


ae) os a Day AW; O55 AW, 





Oates), 
From (11) 
> a. & . 
(S50 Tes ERIM a ahley NAVE ELESN es eeL 
Oso a 42 O42 Aan 
Q a 
pe VN iy a NT (1.18) 
O42 Gaze 
From (13) 
desy Ars 
AQ Ales Ne 
) Os3 As3 
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Substituting AT, from (19) 


Os, die a 
Aa, —4 an nae +4 — AP. (4.19) 


Substituting AP, from (1.17) and Aa, from (1.19) into 
(1.18) and rearranging ; 


Ax =| 45 _ Steg _ %45-4ss Gor) yp _ / S41 Yse\ 4, 


Qas Ass Qao 


4. O46 Sar. 0. bt AX 4 ang a Gar, D3. AW + As7 os maia2 Ay 
ays We O42 Ake One Q42 W 


cs 





OL, 
AX. b APL -b,, Aws bz, Ax +b, Aw, 
+ by, ATy + O,, AW; — Oy, AW ~ Og ATon (1.20) 
From (8) 
AW .- 22 AX _ 928 AP, + 239 AW, 4 <3! AW (1.21) 
Qag Q2g Qag Uag 


o 6 
Substituting AX from (1.20) and AP, from (1.17) into 
(171.21) and rearranging :; 


41 























° a .b Q b. a by Qiee a 
AW, aa |{ 2736 , ~29-"23 AP. 4 27.937 23 eae ai AW 
Gas 229 Gag A293 Qag 
a Jb Qo a. ib a Qa Qa ; F 
= 27° ~38 + 2% bat AX 4 g3° ¥3 36 27 b39 AW 
QA2g A 2g QAog Q29 X29 W 
97 D ho 23 bs 3 AT + CE raloley Qo, by, AW 
= bs a =n l 
Qog Q2g 29g a 29 
sar b + ; ~ ; d b 
ae = 42 _ Gag b 5. AW, + | 27 43 AT 
Usg 123 Qog BW 
OL, 
AWy = - ba, AP, + bas AW - bas AGE De One ATw 
+ bag AW + b, AW, + 6.2, ATgy (1.22) 


Substituting A Wy, from (1.22) into (14) and rearranging 


AP, » (Aen =] AM, 5 (aon ae at} Ont ens i 
| 57 


As57 Os | 


Qs7 


AW; 


ee 


Osa. bss 
BME ae ooo ary 


1 Osg. Dag Ke O55. bag AT, N G5, bag 
Asz 
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Or, 


APy = b, AW, ~6,, AP 4b, AW-b.. AX - bee At, 


+ bs, AW, + bs, AW. 4 b., AT. (1 23) 
From (9) 
AW 2 92 AP 4 239 Ap, _ Se Aw _ 430 Ax _ S97 AW, 
a. . a, . O35 Qys5 Ass 
5 5 


Substituting APy from (1.23) into (1.24) and rearranging ° 








> be 1G Q Cote 
A Ni | SS 2S | aw + (iat — Sts 


35 Q35 A35 Qs, 


OL, 


AW = be, AwWy + be, AP, +b, Aw — Sy ie le 


+ bes AW, + beg Me eye Ne CS) 
From (10) 
: O 44 As9 
AT ay = ae ae AQ. (1.26) 
40 40 ; 
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From (12) 
Qa AB ge = 5 ATI + 85, AWE + O,, AW, - Oy, ATs Cie) 


Since Bye . £rom (21),then substituting AT 





q3 gs 
From (1.9) into (1.27) and rearranging :; 
Aa E Qs, EMiag 2 Qa3.be Ag + ( “322849 a AW 
ye ale ar A4% . 
Q4g.b : 
—(———+ | AT,,, (4-23) 
Cag 


Substituting A ag, from (1.28) and AQ. from (1.19) into 
(1.26) and rearranging : 


. Q..,a Q oa As 
‘ae Speen Se ae Oe Ra. 3 Q.,~ 4 .D. AW 
Q4o-G4g GheGss} PN [a,,-a,, ie u 














| Asg [ass + A4a-b,) AW ~( SesSe 2) an, 
a (3 Gamay © 
4 33: Bee 97 AVES 
Aso: Ass 
Or; 
AT, = Des AT ay . be AW id ee AW, - DS AT si by APs (1.29) 
From (7) 
0, an 
Ave ae eee NP O25 A (1.30) 
9) An 3 X23 A23 8 
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Substituting AP, from (23) and AP. from (20) into (1.30) 


Aw . Qa6 a Qos, Ags AP 2 | AT - G24. aa A 
B An, Aes B Qas ° Ks 


Or, 


th @ Wig MIA = oon Bn See AP. 
From (18) 


AW. = Oe a Ag + sa: AT. 


S 
faa 


Substituting AP. from (23) into (1.32) 


(Ag5.0 
DV areas aia ieee ays SST San. ais dist AP, 
ag5 G93 Ag3 


OL, 
= Oe AXy + Dex AT. + oe APs 


Substituting AW, from (1.31) and AW. from (1.33) 
(1.1) and rearranging : 


APs =(6,.b, by 5) AP, + (b,-b, - oe) a 
+( by bys ) ap, + (d,-b,,) AX, 
Or, 


APs~ Cx BP, + Cy, AT, Cig AP, + DY AX, 


) 


@2ey 9 


(4290) 


(1.39) 


into 


(A) 


Substituting AW, from (1.31) into (1.4) and rearranging ; 
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ae “ by bys APs + ( by + by- bg) AT, - bs ATsy ~ by bya AP, 
QI, 
AT; = Co APs + Caz ATs + Cag ATsy + Cag SP Gey 


Substituting AW, from (1.31) into (1.11) and rearranging ; 
ATs, = ~ is.brs AP, - ( bis. byy + big) ATs + bia ATsw 
Gigs te AVR BINA AWs + by, AW 
OL, 
ATey = Ca; AP, + Cap ATs + Cas ATey + Cag APp 
+D,,44f + D,, Aw, (C) 


Substituting AW, from (1.31) into (1.20) and rearranging ' 


ie oe yp be APS + nor, Bes bs, 4X-b, AW 
+ d,, ATy + D4, AW] 
OL, 


AX 2 Cy, APs + Cay AT, + Cyg AX + Can AWt Cae AWy 


+Ca, NU ee ap GaAs + Cag ATy + Dy, AW; (D) 
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Substituting AW, from (1.31) into (1.25) and rearranging : 


AW = 7 Gee oe APs ia oe : bo. AT. — be, AX+ bo. AW + bs, AWy 


+ he AT aw it ( 5) ya b. + be: | A Ps ~b,, ATw + be AW; 


OT, 
NA a oa mG a Ai 6 pax th Cas AW 4 Csg AWy 
+C57 ATay +Cag AP, + C5q AT, + D54 AW, (&) 


Substituting AW, trom (1.31) into (1.22) and rearranging : 


° 


AW, = ~ 5-6, AP, ~ bg). by, AT, — ba, AX+b,, AW 
+ by, AW, + O51 ATay +( bey bys ~ by) A Py 
- bag ATy + Das AW} 
OT, 


AW, = oy, APs a5 Oe Aare + Coy aS are CEs AW + Cog AWw 


+Coy AT, + C68 A Pe + Cea ATw + Deg AW; 


oo 
“y) 
a 


From (1.29) 


Alpw = - O Maser t oe EN + a AF, + DG A We + DAW, 
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ATan = C73 ATsy + C77 AT, + C75 4 P40, 4M +0,,4N. (GC) 


Substituting AW, from (9.31) into (1.17) and rearranging ° 
AP, 2 ~ bas" b, AP, = bs ba AT. uv oe Ax + Dao AW 
= b,, AW +(b,4 + bys 73) aie t Daa ATy~ 0, AW; 
OL, 
AP. = Cer AP, + Cen ATs + Coy AX4 CQ, AWH ge AWy 
+a, AP, a Cog ATy + Dey AW; CH) 
From (1. 15) 


Me b,, Ax +b, AWw+ bo, 4W, 4b AP +b ATy + bag AW; 


OT, 
ATy = Cg, AX + Cgs AW 4log AW, + Cag AP, +CggATy+0,,AW; (T) 
From (1.13) 
e 
AY = by AX + bg AW + Oa, AW, - Oy, AP, + b,, ATW + bye AW; 


oc, 


AY = Cog AX + Cys AW 4 Cog AWw *Crog APKC, oT, 40 aH, CT) 
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The state variable equations from (A) to (J) were rewritten 
in matrix notation as shown in Figure 4. The state equations 
were soOived for various percent step changes to input 
Variables ( AX, AW, AW,) ,using the IBM simulation language 


csmMP ~M as shown in the next section. 


4g 





XITGLWH ATAVIYVA ALVIS 


0 AV Q 4 
O01} |; LlVii o 
0; avi o 
0; Pivil o o 
0, MVT; @ °9 


7 
| 
| 
: 
| 
| 
| 
! 
| 
4p 
pire 
0 
oOo! AVI} oo 
Oo 
sl 


- 4 Banh ty 


ee) a eee 
ae 0 a 
8 M05 
a os 
as alg Jos 
BS 8) M85 
BUS) ats) a5 
> 0 0 
a (2 0 
ols 0 ” 


Q0F 


by 4 


ras} 














a av 
0 - v 
1) "AV! 
l=] 
re AY | 
ie | XY 
deol | Ly 
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IilI. COMPUTER RESULIS 


ab cm EE Ga ee a cae ee ey ce ee ee 


From Appendix A, it is seen that the a coefficients are 
in terms of steady state values. These values were 
determined from the data in Whalleys thesis (2) which are 
repeated in Appendix E. A Fortran IV program to find the 
a,b,C,D coefficients which appear in the CSHP program is 
shown in Appendix F,and the calculated values are shown in 
Appendix G. Appendix H is the IBM simulation language CSMP 
program that was developed. Only input variables were 
changed for each computer run. For example, results for a 5% 
step change of Ax,, a 10% step change of AWrand a 10% step 
Change of AW; were calculated. The output response curves 
were plotted using the NPS CALCOMP PLOTTER. | 
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A. TRANSIENT RESPONSES FOLLOWING A STEP CHANGE IN THROTTLE 
SETTING OF 5% 


At a particular steady state operating condition the 
throttle valve is suddenly opened while still keeping the 
previous air and fuel flow rates and feedwater flow rate 
constant. Because of bubble formation in the steam drum 
liquid,the steam drum water level (FIGURE 9) swells for 
about 40 seconds and then Shows a steady falling in level at 
higher flow rates of steam. The steam drum pressure (FIGURE 
8) immediately begins to fall. Following these sudden 
increases, the steam mass-flow rates from the steam drum and 
the superheater (FIGURES 10 and 11) show slight declines 
because of the effect of the decrease in st¢am drum 
pressure. The Superheater outlet pressure (FIGURE 10) acts 
in a Similar fashion to the steam drum pressure, that is ,it 
shows a steady decline to a new steady state condition but 
it is slightly lower than the Steam drum pressure because of 
system pressure drops. The superheated Steam temperature and 
Superheater wall temperature (FIGURES 5 and 6) both decline, 
although the latter temperature is slightly lower. The 
riser tube wall temperature (FIGURE 8) and steam drum liquid 
temperature (FIGURE 9) also decrease due te the throttle 
Change. The superheated steam density (FIGURE 5) Shows a 
rising characteristic to counteract the loss in superheater 
pressure and temperature and the gain in steam mass-flow 
rate from tne steam drum (FIGURE 11). The riser and 
downcomer flow rates (FIGURE 7) are reduced to a lower value 
but steam guality (FIGURE 6) is increased to a higher value 


than previously to preserve the energy palance. 
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B. TRANSIENT RESPONSES FOLLOWING A STEP CHANGE OF 10% IN 
FUEL FLOW RATE. 


The effects due to a change in the fuel flow rate are 
less than those due to a change in the throttle valve 
setting. However,the system requires a longer period of time 
to reach a new steady state condition. Because the 
combustion rate is increased, the evaporation rate is also 
increased making both the steam mass-flow rate at the 
superheater (FIGURE 17) and the steam mass-flow rate from 
the steam drum {FIGURE 18) increase in a similar manner. 
Initially, density in the risers falls quickly in response 
to the increased firing rate which causes a rise in 
superheated steam density (FIGURE 12). An increase in steam 
flow from the steam drum without an increase in the 
feedwater flow rate results ina drop in steam drum water 
level (FIGURE 16).The steam drum pressure (FIGURE 15) and 
Superheater outlet pressure (FIGURE 17) rise monotonically 
as time increases. Also superheated steam temperature 
(FIGURE 12) , superheater wall temperature (FIGURE 13) , 
riser tube wall temperature (FIGURE 15) and steam drum 
liquid temperature (FIGURE 16) all increase with increasing 
fuel flow rate. The abrupt fall in riser mixture flow rate 
(FIGURE 14) is probably due to the sudden drop in the level 


of the initial density change in the riser. 
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C. TRANSIENT RESPONSES FOLLOWING A STEP CHANGE OF 10% IN 
FEEDWATER FLOW-RATE 


The main effect of increasing feedwater flow rate while 
Maintaining throttle setting and air-fuel flow rate constant 
is an increase in the steam drum water level (FIGURE23). 
Because fsedwater is added into the steam drun, steam drum 
pressure (FIGURE22) suffers a snall overshoot followed by a 
decline. As a result, a drop in both steam mass-flow rate 
from the steam drum (FIGURE25) and steam mass-flow rate at 
the superheater outlet (FIGURE24) occurs.Flow around the 
riser-downcomer loop (FIGURE21) is only slightly affected. 
Since the firing rate is unaltered,the riser tube wall 
temperatura (FIGURE 22) and the stean drum liquid 
temperatures (FIGUREZ2 3) decline. The superheated steam 
temperature {FIGURE19) and the superheater wall temperature 
(FIGURE20) xcise because of a drop in steam nass-flow rate 
from the steam drum.The quality of steam (FIGURE20) shows a 
general lowering. The superheated steam density (FIGURE19) 
Shows a decline which is similar in form to the drop in the 


Superheater outlet pressure. 
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This analysis is only the initial step toward a better 
understanding of a naval steam generating plant both fron 
the dynamic and from the control points of view. The 
analysis concerns a boiler system which is divided into four 
parts. The major difficulty in the analysis is the fact 
that the whole system is very complex and contains numerous 
variables which are extremely unwieldy to manipulate. Large 
humbers of relationships are non-linear so it is necessary 
to generate a linear form by the methods of perturkation 
theory.After setting up the state variable equations, the 
IBM simulation language CSMP-I was used to solve for the 
open loop transient response characteristics of the state 
variables due to variations in input variables. Throughout 
these calculations,all the variables appearing in the 
transient response curves are not absolute values. 
Rather,they represent a small ircremental change from some 
particular steady state operating condition. Also,it is to 
be noted that trends in some of these results differ 
markedly in form from those appearing in Whalley (2). Future 
studies will be undertaken to determine the causes of these 
differences. Meanwhile,the benefit one may derive from this 
basic analysis is its use in the design and constructicn of 
a new type of boiler controller, that is, a multivariable 


control systen. 
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